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Abstract 

The phase behavior of a melt of monodisperse rod-coil diblocks is studied. We 
derive a Landau free energy functional for both a compositional and a nematic 
order parameter. The excluded volume interaction between the rod blocks is 
modeled by an attractive Maier-Saupe interaction. The incompatibility be- 
tween rod- and coil blocks is modeled by the usual Flory-Huggins interaction. 
For a large volume fraction of the rods, a transition from isotropic to nematic 
to smectic C is observed upon decreasing the temperature, whereas for small 
rod volume fraction, spherical, hexagonal, and lamellar structures prevail. In 
the smectic C phase, the rod orientation angle with respect to the lamellar 
normal increases rapidly from 35 to 40° close to the nematic/smectic-C phase 
boundary to values between 45 and 55°. 

I. INTRODUCTION 

Due to chemical incompatibility flexible AB diblocks are known to phase separate on a 
microscopic scale characterized by the radius of gyration of the blocks. The morphology of 
the microphases in the melt is successfully described in various approaches, which treat the 
AB diblocks as Gaussian chains and the incompatibility by a Flory-Huggins interaction |l]-[Jj . 
Rod-coil diblock copolymers are a special class of AB diblock copolymers. Besides the usual 
incompatibility between distinct blocks, the liquid crystalline behavior of the stiff blocks 
can have a profound effect on the phase behavior. The combination of phase separation and 
orientational ordering gives rise to a wealth of different phases. In this paper a theoretical 
description of both phase separation and orientational ordering is developed. 

In experiments on various rod-coil or liquid crystalline polymeric systems interest- 
ing phase behavior and micro-structures are observed [pp]14l. Linear rod-coil multi-block 
oligomers can have orientationally ordered lamellar, hexagonal, spherical (bcc), and bicon- 
tinuous cubic liquid crystalline structures when the rod volume fraction is less than the coil 
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fraction 0-§J. Schneider et al. reported results on a particular rod-coil diblock copoly- 
mer melt in which microphase separation occurs for higher temperatures than the liquid- 
crystalline order; although at high temperatures the microphase is bcc or gyroid, whenever 
anisotropic orientational order starts to play a role, the lamellar or smectic phase appears. 
One of the most striking features is the appearance of smectic C phases, and so-called zig-zag 
and arrowhead phases [|lOl , |Tl"| for intermediate rod volume fractions. 

The theoretical description of the phase behavior of rod-coil diblock copolymers requires 
the introduction of a compositional (density) order parameter and a (anisotropic) orienta- 
tional order parameter or nematic order parameter. In the work of Holyst and Schick |T5| . 
the interactions and the excluded volume effects such as the steric repulsion between the rod 
blocks are treated in an approximate manner via the Flory-Huggins interaction, Maier-Saupe 
interaction, and the incompressibility constraint. The stiff part of the diblocks can be mod- 
eled by rigid rods il5fl , wormlike chains [|16|] or freely jointed rods |L7||, having a Maier-Saupe 
interaction [I8|. The isotropic-nematic transition in various polymer systems was studied 
in a number of papers [|16| ,|19|-|2T| . Using self-consistent field theory, the interplay between 
phase separation and liquid crystalline order has been studied ignoring morphologies other 
than lamellar f22| , |2~3|| . Gurovich found that copolymer melts can order in four distinct ways 
under the influence of external ordering fields {e.g., electric fields) close to the spinodal point 
[24]. 

For large rod volume fractions, the smectic behavior of layered rod-coil diblock copoly- 
mers in the strong segregation limit was studied p5|- p7| , p3| . For small rod volume fractions, 
so-called hockey puck phases might appear in the strong segregation region [28| or when- 
ever phase separation occurs prior to nematic ordering [29]. Recently, a two-dimensional 
dynamic mean field model using the Maier-Saupe interaction was employed for a melt of 
semiflexible polymers [|21|| . Also the time evolution of morphology formation was studied 
in two dimensions for liquid-crystal/polymer mixtures showing a large influence of nematic 
ordering on phase separation [RO] . Moreover, in two dimensions, the phase behavior of rod- 



coil diblock copolymer melts was obtained using a self-consistent field lattice model [31 



Interestingly, evidence for the existence of metastable zigzag structures which have been 
observed in experiments [|l(],|Il|], was found pif . 



In this paper, along the lines of Holyst and Schick [ 15] and Singh et al. |L7| , a monodis- 
perse incompressible rod-coil diblock copolymer melt is studied with both a Flory-Huggins 
interaction and a Maier-Saupe interaction. The Maier-Saupe interaction models the steric 
repulsion between the rods, and henceforth favors alignment. Special attention is paid for 
the interplay between nematic ordering and microphase separation. One of the main ques- 
tions addressed is: "How does the nematic ordering affect the microphase separation?" The 
Landau free energy is calculated up to the fourth order in two order parameters in the weak 
segregation limit. A density order parameter describes the tendency to microphase separa- 
tion, whereas a nematic order parameter describes anisotropic alignment of rods in the melt. 
For blends, the free energy up to fourth order was already obtained by Liu and Fredrickson 
|nj. Perhaps another way of addressing the current issue is using the density functional 
formalism, e.g., see the paper by Fukuda and Yokoyama |32| and references therein. 

The setup of the present paper is as follows. In the next section a description of the 
partition function describing the rod-coil diblock melt is given, followed by an introduction 
into the Landau mean field expansion. Subsequently, the free energy expansion up to sec- 
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ond order is briefly reviewed, illustrating the instability of the isotropic phase with respect 
to nematic ordering and microphase separation. Then, the actual Landau free energy up 
to the fourth order is derived in the so called first harmonics approximation (FHA). The 
minimization of the free energy and the corresponding phase diagrams are discussed. In 
Appendix [A] the second order Landau coefficient or inverse "scattering matrix" is computed. 
Finally in Appendix |B] and |C| the pertinent three- and four-point single chain correlation 
functions and vertices are given. 

II. THE MODEL 

We consider an incompressible melt of n monodisperse rod-coil diblocks in volume V. 
The flexible (coil) part of the diblock is modeled by a Gaussian chain, and the rigid part 
by a "thin" rod. Each diblock molecule consists of N segments of which Nq = fcN are 
coil segments and N R = f R N are rod segments, with f R = (1 — fc). Thus, the volume 
fraction of the coil part is given by fc = Nc/N. We keep the total density constant, 
Po = nN/V . The characteristic length scale of the coil block is given by the radius of 
gyration, R g = b^jNc/6, with b the statistical coil segment length scale. The length of the 
rod is characterized by i = N R b', with b' the rod segment length. In this paper, we set 
b = b' = 1, thus po = nN/V = 1/b 3 = 1. 

The Hamiltonian for our melt contains two terms: one describing the standard Flory- 
Huggins repulsion between the rod and coil blocks and a second term describing orientational 
ordering of the rod blocks in the melt. The Hamiltonian for our model is 

Hi = x J ' dxp R (x)p c {x) J ' dxS%(x)S%{x), (/i, v = 1, . . . , 3) (1) 

where pc and p R are the densities of the coil and the rod blocks, respectively, and S R is the 
nematic order parameter tensor. The tensor S R is symmetric and traceless. x is t ne usual 
Flory-Huggins parameter. The second term in Eq. ([D) is the Maier-Saupe interaction with 
the parameter uj. It effectively describes the excluded volume interaction between the rod 
blocks in the melt favoring their alignment. 

The conformation of a diblock in the melt is given by the vector function r(r) describing 
the contour of the coil, the vector R giving the position of the joint, and the unit vector u 
describing the orientation of the rod, see Fig. [I]. Then, the (single-chain) partition function 
of a rod-coil diblock in external fields is 

Q[w R ,w c , W^} = C J VfdRdu5(f(f c ) - R)5(\u\ - 1) 

x exp J^ c dT \ff + J WrPr + J Wcpc + J W ^s^. (2) 

where W^ u is an external tensor field which couples to the nematic-order operator. The 
normalization constant C is chosen such that Q[0,0, 0] = 1. The density and orientational 
operators are 

p R (x) = ^ffds6(x-R-su), (3) 
o 
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fc 

p c (x) =N I ds5(x-f{s)), (4) 



o 



Sr(z) = ~g — J S(x — R — su) 



i 



u»u v - 



(5) 



From Eqs. (0)-© all single chain correlation functions can be obtained by differentiating Q 
with respect to the external fields. 

The partition function of the whole incompressible melt can be described by 

2 = | II / ^mdRmdUrnVllf^^Urn}]^^ ~ p R (x) - p c (x)) 

x exp j-x j dxp R (x)p c (x) + ~J dtfSgfflSgfflV (6) 

where 

V[{r m ,R m ,u m }] =C6(f m (f c ) - Rm)6(\u m \ ~ 1 ) ex p{~2^2 J rfr|f m | 2 |. (7) 
After a number of Legendre transformations, we obtain 

Z <x J Vi, R Vi) C VS%VJ R VJ c Vr" 6(i> R + ifc) 



X 



exp {-*/ i'RiPc + ^ J SrSr 
exp {- J J R ^ R - J Jc4>c - J J^Sg - G[J R , J c , J""]}, (8) 



where 



G[J R , Jc, Jn = -n\nQ[J R , J c , J»% (9) 

and where we have introduced the fields 

4>c(x) = p c (x) - f c , ip R (x) = p R (x) - f R . (10) 

Making use of the above definition of G, the Landau mean field free energy can be obtained as 
an expansion in powers of the concentration profiles ipR, ipc and the nematic order parameter 

III. THE LANDAU MEAN-FIELD FREE ENERGY 

We can write the "entropic" and "interaction" parts, respectively, of Eq. (|8]) as follows: 
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exp {-Fj [</>]} = exp {-I a p J 0«0/3 j, 

exp {-F 5 [0]} = JVJ exp |- | J a Q - G[J]}, 



with the shorthand notation, 

p a = (pr,Pc,Sr)- 
Thus the Greek index a = R,C, S. So that 

Z = exp {-F} oc / exp {-F[</>]}, F = + F s 

The entropic part of the Landau free energy reads 

rex 



1 

fel k 2 



+ -^71 11ml <Pa{ k i)<j)f3{k2)<j)^{k3)<j)s{- k i ~h- h) 

k ± k 2 k 3 

x r ( ah 5 (k h k 2 ,k 3 ). 

The vertices Y are 

r§(k 1 ) = [w$j(k 1 )Y 1 , 

Y%(h,k 2 ) = -Y^i^Yf^Y^i-h - k 2 )W% Y (k u k 2 ), 



and 

r2^(fci, k 2 , h) = -r« (fcor^c^r^c^rgc-fc! - a* - fc 3 ) 

- 5 K {h + fc 3 )Wi5(A; 1 )^ ) (A ;2 ) - <fe(A; 2 + k 3 )W® {kjW® (k 2 ) 



- feajrgc-*! - kjw^Uh + k 2 , h 

- W^^h, h)Y^(-h - h)W%,(h + k 3 , k 2 

- wSl^h, -h -k 2 - h)Y$(k 2 + k 3 )W^ lpi (-k 2 - k 3 , h) 
The functions are the single chain correlation functions 

w^ a2 ... an {K h, ■ ■ ■ , fcn-i) = N- n {p ax {k x )p a ,{k 2 ) ■ ■ -p an (- £ h)) 

i=l 

where the average (• • -)o is defined as 

(F[{f,R, u}])o = J VfdRduV[{f,R,u}]F[{f,R,u}], 
where V is given by Eq. (0). 
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IV. THE SPINODAL 



First, we consider the Landau expansion up to second order in the density and nematic 
order parameters, to analyze the instability of the isotropic phase with respect to microphase 
separation and nematic ordering. Therefore, following the approach of Holyst and Schick 



Toj] and Singh et al. [[H]], we consider the free energy 

( m \ 



= l^(-q),S(-q))T^(q) j , (23) 

where r^ 2 ^ is a 2 x 2 "scattering matrix". This matrix can be obtained from the matrix 
r given in Eq. (|A14|) in Appendix |A]; 




r (2) = ( h RR + h cc - 2h RC - 2Nx 2(h RS - h C s)/3 \ (9d) 

[ 2(h RS -h cs )/3 h S s-2Nu/3 J' { ^ } 

where 

see Appendix For the nematic order tensor S R we have taken the Ansatz 

Sg® = S(q) (Xf - ^ , (26) 

where S is a scalar function. The nematic ordering is taken parallel to the wave vector q. 

The "scattering matrix" depends on q, f R , N and the ratio r — u/x- The spinodal 
is determined by the root of the determinant of r^ 2 ^; 

detr( 2 ) = x(q,fR,N,r). (27) 

The real spinodal is given by the global minimum of x (ie., the highest T) with respect 
to q. Whenever one of the eigenvalues (or both of them) of r^ 2 ^ becomes negative, the 
isotropic phase becomes unstable. The onset of the nematically ordered phase is given by 
the minimum of the root x a ^ wave vector q = 0. The onset of microphase separation 
is given by the minimum of the root x a ^ some nonzero q = q*, which gives the inverse 
characteristic length scale of the microphase. The minimum at q = for nematic ordering 
gives rise to a global nematic order parameter S(q) — > S. All this is discussed in detail by 
Singh et al. [JL7|J , although for slightly different microscopic models. The situation where 
spinodal instabilities with respect to two distinct wave vectors occur resembles the situation 
encountered in comb-coil diblock copolymers |33|.|3l . 



V. THE FIRST HARMONICS APPROXIMATION 

In the previous section, we pointed out that the instability of the isotropic phase with 
respect to nematic fluctuations gives rise to a global nematic ordering. The density order 
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parameters are constrained; the volume integral of ipc and i/jr should vanish. Therefore, 
the Fourier components ipc{p) an d ipiiip) can only give contributions for nonzero p. The 
nematic order parameter is not constrained, and it predominantly describes global ordering 
corresponding to the zero mode p = in Fourier space. 

The Landau free energy given in Eq. ( |TTD is not solvable as it is. A way to solve the mini- 
mization condition for Eq. (|17D is to expand the order parameters in harmonics representing 
certain discrete symmetries. Therefore, we shall adopt the first harmonics approximation 
(FHA) for the density order parameter ip = ip^ = —ipc, *- e - ; we take in momentum space 

Hk) * -4= E eM*k{Q - k), \Q\ = q\ (28) 

with q* the radius of the first harmonics sphere, and where E^(Q) is a phase factor and H $ 
is the set of wave vectors Q with radius \Q\ = q* describing a particular discrete symmetry 
S of the structure. The number n$ is half the number of vectors in H$- For the lamellar 
morphology rii am = 1, for hexagonal n hex = 3, for bcc n bcc = 6. Moreover, the phase 
factor E^,(Q) = 1 for the basic morphologies of interest. Clearly, the FHA restricts the 
microphase morphologies to morphologies which can be described by a single amplitude 
or order parameter ip. Especially the nematic ordering could give rise to more complex 
morphologies which require more than one amplitude. To keep the analysis manageable, 
morphologies, which would require more than one amplitude, are not considered in this 
paper. 

As explained in the previous section, we can assume that the nematic ordering is global 
and thus described by a space independent order parameter proportional to the nematic 
director rf'rf — 5^/3. Thus the Ansatz for S^ u is 

S^(k) ~ S5 K (k) (r,^ - 6 ~Pj= S5 K (k)AT^, (29) 

/ \ 2 

= I r/V --§-)' M VV M VV = g- (30) 

where S on right-hand side is the space independent order parameter, and \fj\ = 1. With 
the Ansatze fl28|) and (p9|), the free energy can be expressed as 



eg - 2N X ] ^ + 



' (2) Nu 
C SS 



+ c \h]}ibih' l l )i + cLiibsi' 3 S + c \hl)Ss' l l )2 S 2 + c ssssS 4 , (31) 



where terms of the form ipS, ipS 2 and ipS 3 are absent in this approximation. The c coef- 
ficients are given in terms of 7 functions, which are related to the T functions defined in 
Sec. [Til] (see below): 

i2! - ' EEW + ft)78(Qi), (32) 



2n S Qi VJ 



S sl = ^S(O), (33) 
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cgk = 3,^ -7= E E E 8k(Qi + Q2 + QshfUQi, Q2), (34) 



Qi Q2 Qa 

4ls = ^- s T,J2 8 K (Qi + Qafr&CQi, Q2), (35) 



and 



41 = ^5(0,0), (36) 



C SU = J^Y,EEi: WQi + Q2 + Q3 + Qdl^iQi, Q2, Qs), (37) 

^■ n S Qi Q 2 Q 3 Q 4 

c^s = 1 E E E 8k(Qi + Q2 + Q*h$b s {Qi, Q2, Qs), (38) 

J - nS V n S Ql Q 2 Q 3 

,( 4 ) _ 1 ^ 1 n-^W 

4n s 



4lss = t-EE 8 K (Ql + Q2h^ SS (Ql,Q2,0l (39) 



C^55 =^55(0,0,0). (40) 

In the above formulas the 7 function are related to the T functions of Eq. ( |TTD in the following 
way. The second order 7 functions are 

7 $(Q) = e a e b T$(Q) = T$ R (Q) + T$o(Q) ~ ^(Q), (41) 
7^(0) = M^M pa Tf s ^ pa (<S) = AWg(O), (42) 

with given in Eq. (|30D, and where ipR = e^ip, ipc — ecVs 

e R = 1, e c = -1. (43) 

The third order 7 functions are (with \Qi\ = \Q 2 \ = \Qi + Q2D 

7&(Qi, Q2) = ^^(Qx, Q2), (44) 
7&(Gi, -Qi) = eae^T^iQ,, -Q x ), (45) 
ifssM =NNNTf ss M. (46) 

Finally the fourth order 7 functions are (with |Qx| = IQ2I = \Qz\ — \Q\ + Q2 + Qz\) 

IwwiQu Q2, Qz) = e a e 6 e c e (i r^ ) cd (Qx, Q2, Qs), (47) 

7&s(Qi, Q2, -Qi - Q2) = e a e b e c UY^ cS {Q^ Q 2 , -Q x - Q 2 ), (48) 

7Sss(Qi> 0) = e^AfV^siQu -Q u 0), (49) 

7§5s(0, 0, 0) = MNN. A/T^ 55 (0, 0, 0). (50) 

The coefficients c^ 2 \ c^ 3 \ depend on iV, /r, g*, u>/x, X an d the symmetry-group 
(LAM, HEX, BCC). Moreover, the coefficients and also depend on the angle 9 
between the director rj and the orientation of the lattice of the microphase symmetry. The 
characteristic wave vector q* depends on the architecture, N and 
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VI. MINIMIZATION OF THE FREE ENERGY 



Upon minimization of the free energy Eq. (|3"1~D, with respect to the order parameters if) 
and S and the angle 9, we shall distinguish between seven different phases. These phase are: 

1. The isotropic phase (I): ip = 0, 5 = 

Microphases (M): ip ^ 0, S = 

2. Lamellar (LAM) 

3. Hexagonal (HEX) 

4. BCC (BCC) 

5. Nematic phase (N): if) = 0, S ^ 
Smectic phases (S): if) ^ 0, S ^ 

6. Smectic A: 9 = 

7. Smectic C: < 9 < tt/2 

The smectic phases A and C are defined as lamellar density fluctuating microphases with 
either a nematic director parallel to the lamellar director (A) or with a nonzero angle 9 
between them (C). 

As mentioned previously we assume that temperature and architecture are the two dom- 



inant parameters to play with, and therefore, following the reasoning by Singh et al. [17 



we assume that both \ an d u) scale inversely with the temperature T. This means that for 
a fixed ratio r = uj/x and fixed N, we can draw the usual two-dimensional xN vs / "phase 
diagrams" . 

Clearly, in this fourth order Landau expansion the order parameters if) and S should 
be reasonably small, e.g., if), S <S 1. For the usual microphase separation this is the case 
close to the critical point. However, the isotropic-nematic transition is first order due to 
the presence of the third order term in S in Eq. (|3TD . Therefore, we investigate first the 
"weakness" of the first order isotropic-nematic transition. Assuming that we are in the 
region of parameter space such that ip = 0, we are left with the free energy 



rnem 



(2) Nu 



S 2 + CsssS* + c ssss^^ (51) 



where the coefficients are explicitly given in the appendix. The coefficients are 

(2) 5 (3) 25 (4) 425 

SS ~2f R > Csss ~ 21/1' SSSS ~l96f R 

The free energy (|5T| ) has a first order phase transition to a nematic phase when 



(52) 



c (2) 



Nu 



,( 3 ) 12 f . rJ3) 



-nr =^ c ~ ) ss ~^^\~ll\To2r (53) 

^ssss L ^ssss ) JR 
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where uo c is the binodal value.Q The critical value of S is given by 



,(3) 

-sss 



2c 



(4) 

ssss 



5r ~ 51 



0.27. 



(54) 



This is a typical value for S c for nematic liquids f35fl . Thus for large rod fractions fa, the 
nematic order parameter is close to 1/3. It is well established that the Landau expansion 
for coil-coil AB diblocks is highly weakly first order even far from the critical point in the 
N\-f plane. There, the critical value, ip c , for the order parameter ip is at least one order of 
magnitude less than S c {e.g., ip c ~ 1/30). 

Now, suppose that we are in a regime (coil rich) of parameter space where the nematic 
order parameter is zero (S = 0) and consider the minimization of Eq. (|3l"l) with respect to 
ip. The free energy in this region reads 



J2) 

''ijjlp 



'.2 , J3) J}3 , „(4) ,4 



2N X V + c^V 3 + c 



(55) 



The q* was already (and can still be) determined from the minimum of c?2 with respect to 
\Q\. By definition this minimum is q* . The spinodal is then given by the equation 



J2) 



(56) 



In the FHA the free energy for a lamellar morphology corresponding to q* is known |IJ 
to be given by 



r S 



- 2N X + c 



, -I -I if) 5 



(57) 



where (with \Q\ = q* 



c 



(2)lam 



\ E E ^(gi + Q 2 )7S(Q 

QlGlam Q2&a,m 
1 (4) 



C 



{A)lam 



(2) 



(58) 
(59) 



with 7^ and 7$$^,i explicitly given in the appendix. The hexagonal and bcc morphologies 
have free energies: 



rhex 
■^bcc 



(2)hex _ 2N 1 1 2 , Mhex , 3 U)hex , 4 



(2)6cc _ 2/v 1 , 2 , (3)6cc , 3 , 



J4)6cc ,4 

■'ipi/ji/ji/j r 6 ' 



with 



(60) 
(61) 



„(2)/i«s 



(2) 



^(3)hex 



1 (3) 



(A)hex 



„( 4 ) 



(1) 



12 



2" 



(62) 



x The value w s = 15/(2/1) is the so-called spinodal value corresponding to the change of sign of 
the second order term. 
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respectively, 



c 



J2)bcc 

(4)bcc 
ipipipip 



(2) 

7i/n/n 



„(3)6oc 

-'ijjlplp 



3VQ 



(3) 



,(4) 



,( 4 ) 



(4) 



„( 4 ) 



24 



6 



(63) 



For fixed iV"x, r, and AT, it is then straightforward to find the morphology with lowest 
free energy. 

The next step is to consider the region of the phase diagram where ip <C S < 1. In this 
case if the temperature is lowered first the isotropic to nematic transition is encountered 
meaning that the lowering of the free energy is primarily driven by the appearance of a 
nonzero value for S. When the temperature is decreased (N\ T)? the appearance of a 
microphase structure (ip ^ 0) lowers the free energy further. In this particular case, the 
value of S is mainly determined by Eq. (|51~1) . The value of ip then follows from the free 
energy part of Eq. fl3~T|) , 



+ 4l s S + c^ ss S 2 - 2N X ] r + 



,(3) 



,(D 



■0 3 + c^ip^ 4 1 



(64) 



in which S should be considered as an external field. The total free energy being Ts = 
jrnem _|_ jrmix _ Tn ^ e above Eq. fl64|) the angle 6 comes into play. 

For instance, if the morphology is assumed to be lamellar, the coupling term c^ s can 
be written as 



c 



(3)lam 



(3) 
7^5 



2n L 



E p 2(Q-v) 



(3) 
lipipS 



Qelarn 
,(3) 



P2{Q-v) + P2(-Q-v) 



P 2 (cos^) 7l (3) 



(65) 



with Q = Q/q*, cos 9 — Q ■ r\. In case of a hexagonal morphology (with six lattice vectors 
Qi), we get 



c 



(3)hex 



(3) 
lipipS 



2n h , 



E p 2(Q-v) 



(3) 3 



ea; QGhex 



i=l 



(66) 



and for the bcc morphology 



(3) 

4T* C = E w ■ »?) 

Q^bcc 



^^bcc 



(3) 

6 



E^(Qi-i7), 



The other coupling terms in Eq. 



i=i 

are decomposed as 



(67) 



4n 



s V m Qesym 
X 4 ) 



7^SS0 + 7^S51^1s»/m(^) + 7^S52^2sj/m(^) 



(4) 



(68) 



where sym stands for either one of three basic morphologies, and 
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A lsym {V) = - £ P 2 (Q- V ), A 2sym {Q) = - £ P^Q-v), (69) 

An sym Qtzgym ^ n sym Q &syrn 

with the three dimensional space angle between the nematic vector fj and the orientation 
vector of the lattice corresponding to the space symmetry group sym. 
Subsequently, the free energy is written as 

Tf* = ai ip 2 - a 2 ip 3 + (70) 

with 

«i = + tfisS + c^l ss S 2 - 2iV X , (71) 

a 2 — — c\p\^ + c \p\^S^ i ( 7 ^) 

Then for the lamellar morphology (thus smectic), these coefficients can be expressed as 
a x = -2N X + {3o(q*)+f3i(q*)P2(cos6)+[3 2 (q*)P*(cos6) 1 



1 



a 2 = 0, a 3 = ^7v$Wi> ( 74 ) 
where cos# = (j] ■ Q)/q*. The functions /3j are 

Now the minimum of JF™ X with respect to corresponds to the minimum of a\, since 03 
does not depend on 6. Since /?i < and f3 2 > 0, the minimum of a>i with respect to cos 6* is 

ai = -2N X + Po-jl, (76) 

with 



min = arccos ^^jfj ' ( 77 ) 

Numerically it turns out that ^lgg 2 is roughly one order of magnitude larger than ^Iggi 

and 7^550 f° r a l ar g e range of volume fractions fa- Since the nematic order parameter S 
is close to 1/3 or even larger (e.g., see Eq. (0)), the consequence is that in general j3 2 3> (3\ 
and therefore the angle 9 is close to 



9 m « arccos ^1/3 « 54°. (78) 

Thus a smectic C phase with an angle of about 54° is obtained in the region of parameter 
space where nematic instabilities precede density fluctuations. 

However, the above conclusion is premature unless the other two morphologies are con- 
sidered. For these morphologies, besides the terms already given in Eqs. flBTf ), (|67|), and 
(E8|), third order density terms are present and consequently the coupling term 
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c 



(4) 

■ipipipSsym 



J2 £ itwsiQuQii-Qi-Q*)* 



3\n 



syui Q^^sym Q2^sym 



(79) 



plays a role. This term is written as 



J4) 



C 



(4) 

ipipipSbcc 



3V3 
4 

3^/6 



(80) 
(81) 



with Q 3 = —Qi — Qi and o$^$ defined in Eq. QB78|) . The three vectors Qi can be an 
arbitrary triangle in the vector space of the first harmonic sphere of the particular lattice 
symmetry group (HEX or BCC). However, by comparing numerically the free energies of 
the hexagonal and bcc morphology with nonzero S in the FHA, it turns out that the smectic 
C phase (thus lamellar) has the lowest free energy in this particular region of the parameter 
space. 

In the coil rich region of the parameter space, it is possible to take a value for r = uj/x 
so that microphase separation occurs prior to nematic ordering, i.e., with S <C ip < 1, on 
decreasing the temperature. In this case the order parameter ip is mainly determined by 
Eq. fl55|) . However, now the nematic order parameter S is driven by the field ip which acts 
as an external field in the free energy part 



T 



(3) 



i> 2 + c^sip 



s + 



J2) 



+ c, 



(4) 

ipipSS 



4> 2 - 



(82) 



up to second order in S. We assume that the temperature is still above the binodal temper- 
ature (thus below Nx c ) so that nematic ordering is purely density driven. Then 





r (3) J,2 




2 


\c {2) 
c ss 


- Nu/3 





(83) 



For a large region of the phase diagram with the above constraint S <C ip < 1 it turns out 
that S m <^ S c as given by Eq. (j54|) . This concludes the analysis of the free energy, Eq. (|3lD . 



VII. THE PHASE DIAGRAM 

In the previous section, we outlined the various steps in the minimization of the FHA 
free energy Eq. (|3l|) . This section is devoted to the actual computation and derivation of 
the phase diagrams. The first step in the process is to determine q*. This wave vector q* 
is given by the minimum of 7li(<z) of Eq. (|B9|) with respect to q. It is to be expected that 
the characteristic length scale for microphase separation is predominantly determined by 
the length I = Nfn of the rod part of the diblock. At least for large N and in the rod rich 
region this is to be expected, since the length scale of the rod part scales with N and that 
of the coil part with A 1 / 2 . This means that the lowest q vector for the rod is much smaller 
than the lowest q vector of the coil part. The corresponding characteristic wave vector for 
the rod is q* = 2n/£ = 2tt /(Nfn). In Fig. |] the numerical value of q* is depicted vs rod 



13 



fraction fa and compared with the above-mentioned characteristic wave vector. From this 
figure it is clear that for a large region of both q* lie rather close together, supporting 
the view that the rod-length scale characterizes the microphases. 

With q* numerically computed as function of the higher point vertices and 7^ 
are determined, since these do not depend on x and to. Subsequently all c coefficients of the 
Landau free energy Eq. (|3l|) are computed. Depending on the values of Xi r ; N the phase 
diagram is constructed by minimizing the free energy with respect to ?/>, S, and 9. 

As pointed out in the previous section, when there is no Maier-Saupe interaction, there 
will be no accountable or notable value for the nematic order parameter S {i.e., S <C 1/3). 
The phase diagram for this situation where uj = 0, i.e., r = 0, is given in Figure |3| for 
N = 40. The asymmetry of the phase diagram with respect to the AB coil-coil diblock 

phase diagram (as first derived by Leibler M) is apparent. First of all, the critical point, 

(3) 

which is given by the root of the 7^^ vertex, is shifted to the rod rich part. Beyond the 
critical point for even higher rod fraction the bcc phase is absent. 

The coil rich part of the phase diagram is quite similar to an AB coil-coil diblock phase 
diagram with first the appearance of a bcc phase prior to hexagonal and lamellar phases. 
However, close to fraction fn ~ 0.4 the bcc phase disappears at the triple point and for 
increasing x only the hexagonal and lamellar phases are encountered. This is a result of the 
relative smallness of the ratio k, 



(3)6cc \ 2 (4)hex 

k = I (84) 

1 (3)hex J (4)oec ' 



as compared to coil-coil diblock copolymers [[J. Whenever k < 1 the hexagonal phase is 
encountered first upon crossing the binodal curve from the isotropic phase. The triple point 
corresponds to k = 1. 

The dependence on the ratio r is depicted in a series of figures, Fig. 0-Fig. |7|. The 
isotropic-nematic phase boundary (PB) is roughly given by the spinodal curve 

Xs = 15/(2r/g), (85) 



see footnote in Sec. [VT|. Thus by increasing r, Xs is lowered which appears as a shifting of 
the isotropic-nematic PB to the left. Whenever the nematic phase region "overlaps" with 
the microphase region a smectic C phase (lamellar) is favored over the bcc and hexagonal 
morphologies; nematic ordering favors lamellae. This results from the fact that, whenever a 
nonzero value for S appears, the second order coefficient a.\ in the Landau free energy for 
ip (Eq. (|70l) ) is considerable smaller for the lamellar phase than for the hexagonal phase, 
after minimization with respect to angles. Consequently, in the nematic phase region, the 
spinodal curve for the lamellar morphology (=smectic) lies below the binodal curves for 
the hexagonal and bcc morphologies. In the FHA the hexagonal and bcc morphologies are 
incompatible with nematic order, and a smectic (layered liquid crystalline) phase prevails. 

Contour lines of the angle 9 are depicted also in the figures. The dots denote "triple 
points". The angle increases rapidly from 35 to 40° close to the nematic smectic-C phase 
boundary to values between 45 and 55°. This behavior is in agreement with the analysis in 
the previous section. In Fig. the nematic phase region has completely overwhelmed the 
microphase region. Even in the coil rich region there is an isotropic-nematic transition. This 
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means that the microphases are absent all together and only the smectic C phase appears 
for lower temperatures. 

If the overall length of the diblock is increased (N f) the critical point and the triple 
points will shift to the rod- rich part. However, qualitatively, the phase diagrams for fixed r 
but with different N are quite similar, as can be seen by comparing Figs. (|5|) and (|8]). 

VIII. CONCLUDING REMARKS 

In this paper, we studied an incompressible melt of rod-coil diblock copolymers. Both 
microphase separation and orientational ordering were investigated. We have derived a Lan- 
dau expansion of the free energy of the melt up to fourth order in the two order parameters 
ip and S, representing, respectively, compositional and nematic ordering. The composi- 
tional ordering or microphase separation was driven by the usual Flory-Huggins interaction, 
whereas the orientational ordering was driven by a Maier-Saupe interaction. Up to seven 
different phases of the melt could be distinguished as function of molecular architecture, 
temperature, and relative strength of the molecular interactions. The characteristic length 
scale for the microphases turned out to be roughly the length of the rod part of the diblock. 

In comparison with the coil-coil (AB) diblocks, the phase diagram for rod-coil diblocks 
is quite asymmetric, with a critical point lying in the rod-rich region on the spinodal curve. 
Nevertheless, the spinodal curve for the rod-coil system is nearly symmetric. There is a 
complete absence of the spherical or bcc microphase in the rod-rich region, beyond the crit- 
ical point. This is perhaps not so surprising, since it is hard to imagine spherical micelles of 
coil segments being embedded in a matrix of rods. The coil-rich region is more similar to a 
diblock, where the spherical microphase appears prior to hexagonal and lamellar structures 
when the temperature decreases. However, for intermediate rod-volume fractions, there is 
an interesting suppression of the bcc morphologies with respect to the hexagonal one. When 
crossing the binodal curve from the isotropic phase, the hexagonal phase appears instead of 
the bcc morphology. There is an isotropic-hexagonal-lamellar transition. Such a transition is 
rather uncommon for coil-coil (AB) diblocks. The general features of our phase diagram are 
in qualitative agreement with experiments on rod-coil systems with small and intermediate 
rod volume fractions Moreover, in these experimental papers, observations of bicon- 

tinuous cubic phases are reported. The description of these more complex phases is beyond 
the scope of the present paper, since it would require the introduction of an additional order 
parameter, namely, the amplitude of the second harmonics. 

Most of the features mentioned above are governed by compositional ordering, which is 
most prominent in the coil-rich phase. However, different phase behavior is obtained when 
the orientational or nematic ordering is the driving force in the melt. Therefore, in the 
rod-region, with reasonable values for the Maier-Saupe interactions the main results are; 
the suppression of hexagonal and bcc structures in the rod-rich region of the phase diagram 
whenever nematic ordering occurs and the appearance of a smectic C phase in the rod-rich 
region, with a characteristic angle 45° < 9 < 55°. 

We have shown that when nematic phases appear prior to microphases upon decreasing 
the temperature, instead of the usual bcc or hexagonal morphologies a smectic C (lamellar) 
phase is obtained. Also we have found no evidence for smectic A phases in the rod rich 
region of the phase diagram. As was mentioned in the introduction, smectic C phases and 
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equivalent phases have been reported for experiments on rod-coil systems for intermediate to 



high rod volume fractions ]10||TT|. Moreover, the smectic angle 9 is in quantitative agreement 
with the experimental observed values for 9, with 9 ~ 45°. 

At a first glance, the absence of a smectic A phase in the rod-rich region might seem to 
contradict earlier theoretical results p3| , f25| -|2T|j , where the orientational alignment is chosen 
to be perfect (all rod segments are parallel, i.e., S ~ 1). However, one should realize that 
the present model describes a different limiting case of the rod-coil diblock copolymer melt. 
Clearly the free-energy expansion in order parameters breaks down, whenever one of the 
order parameters becomes too large, e.g., when 5* approaches unity. Roughly speaking, 
the applicability of the Landau free energy is limited to the region in the %A ~ fa plane, 
where the spinodal curves for order parameters ip and S intersect or lie reasonable close 
to one another. Although, it can be shown that the Maier-Saupe interaction is a valid 
approximation for the steric repulsion even when the alignment of rod-segments is nearly 
or completely parallel |T^| , the present model cannot really be compared to the other works 



p3| , |25| -f27| in great detail, due to the breakdown of the Landau expansion for the order 



parameter S. 
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APPENDIX A: THE SECOND ORDER VERTEX 



In this appendix, we compute the second order vertex r^ 2 ^. The vertex r^ 2 ^ is the inverse 
of the single chain correlation function as given by Eq. fll8|), 



(Al) 



where the Greek index a = R,C, S. We introduce the matrix 



( 



(A2) 



where 



Wj£(p) = N- 2 (p R (p)p R (-p)} = 
= f R K$(y) = g RR (p), 
with f R — 1 — fc, and where the variable y is 



'h 



JO 



dsds' (e 1 



(A3) 
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y = N RP . (A4) 

The other functions are 

W&(p) = f 2 c K$(x) = g cc (p), (A5) 
W$(p) = W§l{p) = f c f R K%\x)K R l) {y) = gRC (p), (A6) 

where the variable x is 

x = N c p 2 /6, (A7) 

and Kq\ = fo is the Debye function and the K functions are defined in Appendix y. 
Furthermore, the two point single chain functions with one nematic tensor is 

W R T( P ) = WiT(p) =(ff- 5 -Y) f^nkv) = ^9ns(p), (A8) 

Wgt{p) = Wftip) = - x) MoKS\x)K^(y) = A^gasip), (A9) 

with A Mi/ = p^p v jp 2 — 5 pu /3. Finally, the two point function with two nematic tensors is 
expressed as 

wft pa {p) = fijiK^^rr^p) =j:%^(p)gss l (p), (aio) 

1=1 i=l 

where gssi = f R Ksi \v)i with Kg} defined in Appendix |C[ The tensors are 

T pvp(J (p) = 8^8^ + 8 va 8 w - ^-5 pu 5 pa 

_ 26 v P P^_ _ 26 uaP^_ _ 2 §M pUpa _ 28 mVPpV 

p2 p2 p2 p2 



+ Ur£L + hr&L, ( ah) 

3 3 ;r 



Tr pa (p) = 8 vpP -4r + 8 V ° P -A- + 8» pp -4- + 8 paP 4- 

r p *jp *jp r p 



4 n p n cr 4 nPn u 4 

- ;T - o^jT + ( A12 ) 

3 3 p l 9 





8^\ 


( p p p a 


8 pa \ 


V p 2 


~~3~J 


V p 2 





Tr-(p)^\^r--\[^r--\, (A13) 



Subsequently, the inverse of the matrix W is written as 



/ r (2) r (2) r (2)pa \ 

1 RR 1 RC 1 RS \ 

r pupa _ F (2) F (2) F (2)pcr 

- 1 — L CR L CC L CS 

F (2)/ii/ F (2)/tf> -p(2)fiup(T I 

\ L SR L SC L SS I 



(A14) 



and the following notation is introduced: 
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Wi a) = g ab , r$ = h ab , (A15) 

w m,u = AflVg ^ r (2w = Afa , h ^ (A16) 

= rr^gssn r^- = j2 rr^hssi. (ai?) 

i=i i=i 

The h functions are related to the g functions by 

2 

hab = gab ~ ^9ac9cshs, (A18) 

^5 = -7 J ^ 7, (A19) 

2 (S'cSSW ^dS + $9 SSI - 4^552 - 9SS3 ) 

^51 = ^551, (A20) 

, _ 9sSl9SS2 rA<?n 

*{gss2 - 9SS1) 
hss3 — 7: 3g c sh cS + 32g SS ih S si — ^6gss2hssi 



(A22) 



2 

10(?S53^S51 — 16gsSlhsS2 + &gsS2hsS2 + 8(7553/1552 
1 



5(7551 — 4(^552 — gsS3 

Hence, we have obtained T^ 2 ^. 



APPENDIX B: VERTICES 

In this appendix, all vertices up to fourth order which appear in the FHA are given. 
First, the so called pure nematic vertices are discussed. 

1. Nematic vertices 

The pure nematic vertices are those involving correlation functions of the nematic order 
parameter . Clearly, these are 7^5, 7555, and ^fsss- For instance, the vertex 7 s 2 j can 
be obtained from Y < gg ivpcJ as defined in the previous appendix. 

In what follows, the limit p — > is considered, with p parallel to the orientational vector 
7/. 2 Consequently -> . Thus 



7 g>(0) = M^M p(7 Tf s pvpa (<d) 



A 20 Ifi 4 

^hssiX^M^M^ = — -WO) + -^Wo) + -h SS3 (0) 
,-=i y y y 



2 One can show that parallel alignment of p and ff has highest instability temperature flT^ . In 
other words close to the spinodal, the angle 9 between p and ff is zero. 
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= — lim 



(9cs(p)9cd(p)9ds(p) + ^9ssi(p) - ^9ss2(p) - gsssip)) Jr 



2 ' 



The three-point nematic vertex is 



7§s(0,0) =^^(0,0) = -^(0)] s W^ fl (0,0) 
n 3 

ATAfATW^siO, 0) 



^(3) 



,(2) ( 



r(3) 



15 



15 
2/Ij 



1 3 



/R U 1 35 f R 7 ' 



where we have used the identities 
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p->0f 



i=i 



2/1 



and 



M^N^N^Wf^^^Q) = /!(2/3) 3 ([P 2 (cos#)] 3 ) 

i 

= ^(2/3) 3 / dx[P 2 (x)f/2 = / 3 (2/3) 3 (2/35). 
-l 

Subsequently, the four-point nematic vertex function can be computed, 
t£Ls(0, 0, 0) = AfAfAfAfrf sss (0, 0, 0) 

= -[ATrg(o)] 4 ^S S5 (o,o,o) 

+ 3[ATr^(o)] 4 ^(o, o)ri(o)^ s (o, o) 

+ 3[Arrg(o)]M 2 j(o)^| 2 i(o) 

1 2550 



49 ' 



where we have used that 



AfAfAf, AfW^% s (0, 0, 0) = / 4 (2/3) 4 ([P 2 (cos£)] 4 > = / 4 (2/3) 4 (3/35), 
AfAfW^ s (0,0) = / 3 (2/3) 2 (2/35)A/\ 
A/Wi(0) = / 2 (2/3) 2 (l/5). 



2. Density coefficients 



(2) 



by 



Now the expressions for the pure density vertices are listed. The coefficient 7^ 



7$(Q) = e e 6 r^(Q) = ^(Q) + /^(Q) - 2/ lJK7 (Q), 



,(2) 
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where the Roman indices a, b sum over R, C. At this point, we introduce the shorthand 
notation 

l ab = T^{Q). (BIO) 

The coefficient 7^^ is 

lH,i,{QuQ2) = e a e b e c T ( ^ c (Q 1 ,Q 2 ) 

= -z a (Q 1 )z b (Q 2 )z c (-Q l -Q 2 )W®(Q 1 ,Q 2 ) (Bll) 

with \Qx\ = \Q 2 \ = \Qi + Q2I (thus Qi ■ Q 2 = -1/2) and where 

z a (Q) = e^rgiCg), (B12) 

hence 

zr{Q) = h RR {Q) - h RC (Q), z c {Q) = h RC (Q) - h cc (Q). (B13) 
The three point single chain correlation functions are (again with Qi ■ Q 2 = — 1/2) 

W R % R (Qu Q 2 ) = w RRR , w RRR = f R K$(y, -1/2), (B14) 
W R i c (QuQ 2 ) = w RRC , w RRC = f%fcKg\x)K$(y, -1/2), (B15) 



= WRRR, 


WRRR 


= Wrrc, 


Wrrc 


= wrcc, 


Wrcc 


= wccc, 


wccc 



Wjg c (Qi, Q 2 ) = w RCC , WR CC = f R fhKgl(x, -l/2)K%\y), (B16) 
W^ c (Q ll Q 2 )=wccc 1 w ccc = f c K§}(x, -1/2), (B17) 



where w RRC = w RC r = w C rr and w RC c = w C cr = w C rc- Again the K function are given 
in Appendix |C[ 

The coefficient 7$$^, is (with \Q X \ = \Q 2 \ = \Q 3 \ = \Q t + Q 2 + Q 3 \) 

7$,w(Q 1 > Q 2 > = ^b^dT^cdiQu Q2, Q3) 

= Mi(Qi, Q 2 , Q 3 ) + M 2 (Q U Q 2 , Q 3 ) (B18) 

where 

M 1 (Q 1 , Q 2 , Q 3 ) = -^(Qi)^(Q2)^(Q 3 )^(-Qi -Q2- Q 3 )W^L(Qi, Q2, Qs), (B19) 
M 2 {Qt, Q 2 , Q 3 ) = z a (Q l )z b (Q 2 )z c (Q 3 )z d (-Q l -Q 2 - Q 3 ) 

x wSl(Qi,Q2)r ( 3(-Qi - Q2)W®(Q 1 + Q 2 ,Q 3 ) 
+ w^2(0i.0s)rg(-0i - Q s )wg(Qi + g 3 ,g 2 ) 

+ WSL{Q U -Qi -Q 2 - Qz)vf)(Q 2 + Q 3 )W^ b (-Q 2 - Q 3 , Q 3 ) 

+ 8 K {Q X + Q2)W%\q 1 )W%\q 3 ) + 8 K {Qi + Q3)W^(Q 1 )WS\Q 2 ) 



+ 5 K (Q 2 + Q 3 )W£>(Q 1 )W£\Q 2 ) 

Parts of the form 



(B20) 
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z a z h z c z d W^ s Tf s W^ (B21) 
though formerly present in Eq. (BO), are not included in Eq. (|B2(J|), since these parts can 



be shown to be negligible with respect to the other terms comprising the four point density 
vertex. 

We can compute the four-point single chain correlation functions. Their contribution is 
decomposed as 

W R i RR (Q 1 ,Q 2 ,Q 3 ) = f R K R A \y,c 1 ,c 2 ), y = £q*, (B22) 
W^ RR (QuQ 2 ,Q 3 ) = fcf R Kg\x)K£}(y, Cl ,c 2 ), (B23) 
W& RR (Qu Q 2 , Q 3 ) = f 2 c f 2 R Kgl(x, Cl )K^(y, c x ), (B24) 
wSj% r {Q u Q 2 ,Q3) = fy R KP 2 {x, Cl ,c 2 )K R l \y), (B25) 
W® cc (Qi,Q 2 ,Q 3 ) = f^\x,c u c 2 ), (B26) 



where 



_ Qi ■ Q2 _ Qi ■ Q3 _ Qi ■ Q3 1 /do m 

Cl = loF> C2 = 1qF' C3 = ^T = 1_C2 " ( 7) 



Also, we have (with c = Q\ • Q 2 /\Q\ 



2\ 



W { R 3 UQu Q 2 ) = ftKg>(y, c), (B28) 

W R 3 RC (Q h Q 2 ) = f R f C K£\x(2 + 2c))K$( y ,c), (B29) 

wgUQnQi) = f 2 R f c K£\x)K%(y,c), (B30) 

W^ R (Q h Q 2 ) = ftfcK§\x)Kg>(y,c), (B31) 

W R %(Q h Q 2 ) = M'K^c^iy), (B32) 

W§ RC {Q U Q 2 ) = fR%K%(x,c)K$;\y), (B33) 

^ccr(Qu Q2) = fRfcK { cl(x, c)K R l \yj2 + Yc), (B34) 

W^ C (Q U Q 2 ) = foK^c). (B35) 

For the three morphologies it is sufficient to consider four particular configurations of 
"four vectors adding up to zero" for T^L^Qi, $2; $3)- These configurations are, using 
Eqs. flB2?D, 

ci = -1, c 2 = -1, c 3 = 1 — ► tS^i, (B36) 

d = -l, c 2 = -l/2, c 3 = l/2 — 7^2, (B37) 

ci = -1, c 2 = 0, c 3 = — ► 7^3. ( B38 ) 

Ci = -l/2, c 2 = -l/2, c 3 = — 7^4, (B39) 

where = 7v$w(Qi> W 
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3. Mixed vertices or coupling terms 



The most interesting vertices are those involving both nematic as well as density order 
parameters, since these will give rise to coupling between nematic ordering and microphase 
ordering. 

(3) 

The coupling vertex 7^5 is 

7&(Q, ~Q) = e a e b NT%{Q, -Q) 

= -z a (Q)z b (Q)AfT { p s (0)w£l(Q,-Q) 

= -z a {Q)z b {Q)Nrf s {Q)NN{Z/2)W^ s {Q, -Q) 



-^a(Q)z b (Q)z s W%l(Q,-Q), 



where we have defined 



z s M^M px = r 



(2)fiupX 
SS 



(0) 



45 



zs 



2 • 



and 



W^(Q,-Q)=ATW^(Q,-Q). 
It is straightforward to show that 



W R %f(Q,-Q) 



(Q), 



fnWgt(Q), 



w®TO,-Q) = o. 

wrrsP2(Q ■ v), wrrs 

wrcsP2(Q ■ rj), w Rcs 
0, wees = 0, 

(4) 



Hence 

W R %(Q,-Q) 
W R %(Q,-Q) 

where w RCS = w CRS - 



The four-point vertex 7^55 reads 



(B40) 



(B41) 



(B42) 



(B43) 
(B44) 
(B45) 



(B46) 

(B47) 
(B48) 



7$ss(Q, ~Q, 0) = e a e b ^T$ ss (Q, -Q, 0) 

= - (l)\(Q)*(Q)4[w!&s(Q, -Q,o) - wS(o)wi 2) (Q) 
- -Q)wSs(o, 0) - 2w^ c (q, o)r2 ) (Q)wffi (0. -Q) 



(B49) 



with 
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W$ SS (Q, -Q, 0) = MAfW^ ss (Q, -Q, 0), (B50) 
W^(0)=AfAfW^(0), (B51) 
W^liQ, -Q) = MW^siQ, -Q). (B52) 



For instance, W R 4 RSS is 



and 



Thus 



Also 



and, 



W^ SS (Q, -Q) = flw£t Pa {QW LU M pu = fkE K shvWm(Q), (B53) 



i=i 



rr pX {Q)M^U pX = -yP 2 (Q-v)~ ~, (B54) 
rr P \QW u Af pX = (Q-v)+ |, (B55) 
Tr P \Q)M^M pX = ^ [P 2 (g • V )} 2 . (B56) 



w4rss(Q, -Q) = Wflflsso + wrrss^Q • v) + wrrsszPHQ ■ V), (B57) 
Whhsso = /£ (§)' {-KfKv) + 2K%>(yj) , (B58) 

^{l) 2 (- AK sl(y) + 2Kg(y)), (B59) 



WRRSS1 

^552 = /U^) 2 4 2 3 ) (y). (B60) 



i=i 

= ^ijcsso + w^siPalQ • V) + v>rcss2P£(Q ■ ri), (B61) 
w«75S0 = fcf R Kg\x) Q 2 (-4¥(y) + 21fg(y)) , (B62) 



^cssi = fcf R K^{x) y (-4tf#(y) + 2#&(y)J , (B63) 

P 3 ^(1)/^ f 2 \' z^(l). 



«>*7flW = fcf R K'c\x) [-) K£(y). (B64) 



W& SS (Q, -Q) = flfcK^lix)- = wccsso- (B65) 



The three-point single chain correlation parts of Eq. ([B49) can be obtained from Eqs. ( |B46|) - 
( |B48D and Eq. 
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Subsequently the 7^55 vertex is decomposed in the following way, 



7$ss(Q, ~Q, 0) = 7&SD + itlss^i (Q-v)+ ltlss2 Pi (Q-v), 



where 



(4) 

l^sso — 



(4) 



(4) 



ZrWrrssO + 2zrZ C Wr C SS0 + Z c W C CSS0 



W S s{ZrWrr + 1zrZ C Wrc + Z c W C c) 

(J) Z S Z R W RRSS1 + 2z r z c wrcssi 
zs(z r wrrs + 2z r z c wrcs)w S ss 



ZrWrrsS2 + 2zrZ C WrcsS2 



-2z a Z b W acS lcdWdbS 



Finally, we have the four-point vertex function 7^^ s 



(4) 



(Qi, Q2, Q 3 ) = e^e^itU^i, Q 2 , <5s) 



-g^o(Ql)^(Q2)^c(Q3)^S 

W^ e (Q 1 ,Q 2 )T^(Q 3 )wj%( 
W^(Q 1 ,Q 3 )r^(Q 2 )wjH(-Q 2 ,Q 2 ) 



x 



wSi(Qi,o)rff(-Q 1 )w)S(Qi,Os) 



^(2) 



(3), 



with Q 3 = -Qx - Q 2 . In the FHA, the terms W^W^ vanish. 



It can be shown that the correlation function W <y4; 



rrrs 1S 



W R i RS (Qi, Q2, Qs) = wrrrs [p 2 (Qi ■ 77) + P 2 (Q 2 • v) + P2(Qs ■ v) 
wrrrs = ^f R K { s\y), 



and 



W^RsiQu Q 2 , Qs) + W^rs(Qi, Q2, Qs) + WZcsiQi, Q2, Qs) = w R r C s 



K4) 



(4) 



X 

2 



P 2 (Q 1 -v) + P2(Q2-v) + P2(Qs-v) 



wrrcs = lflfcK^\x)K^ s {y). 



-(2), 



The correlation function W C A CRS is (with Q 3 = —Qi — Q 2 ) 
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W^ RS (Q h Q 2 , Qs) + W^ CS (Q U Q 2 , Q 3 ) + Wg> cs (Q u Q 2 , Q 3 ) = w CC rs 



r(4) 



(4) 



wccrs =~-flflKP 2 {x)K§{y). 
The correlation function Wq CCS vanishes in the FHA. 



(B75) 
(B76) 



We can derive that, using Eq. ( B71|) , ( B73|) and ( B75|) , the four-point vertex (which is 
fully symmetric in its arguments) is 



7^s(Qi> Q2, Qa) = \P2<kQ\ ■ v) + p 2(Q2 • v) + p 2(Qs ■ v) 



(B77) 



with Q 3 = —Qi — Q 2 , and where, 



iptptpS 



2 r 



3 2 2 

ZrWRRRS + z rZcWRRCS + ZrZ c W C CRS 



+ o Z S 1^ Z a Z b Z c W abd JdeWecS, 



(B78) 



a,b,c,d,e 



and zr = z R (q*), z c = z c (q*)- Note that indices a, b, c,d,e — R, C. 



APPENDIX C: DEFINITION OF COIL AND ROD FUNCTIONS 

In this appendix, we list the so called coil and rod functions which appear in single chain 
correlation functions and vertices as given in the previous appendices. The superscript 
of such a K denotes the number of integrations or internal points involved. For most K 
functions an explicit expression can be given, however for certain rod functions the integral 
form is retained, since no analytical form could be obtained. 

The coil functions are 



and 



K§\x) 



1 

x 



/d{x) = f 



e ~ +x — 1 



K<g{{x), 



(CI) 

(C2) 
(C3) 



K<H{x,c) 



Kg}(x) 
Kgl(x) 



1 1 

J dnj dT 2 e- x{1+c)i2 - T1 ~ T2)+xclT1 - T21 



e -2(i+ c )x + (1 + 2c) - 2(1 + c)e~ x 



Kc2^ x i 1 



1 + 3c + 2c 2 )x 2 



-1/2) = 3 



x- 



1 + x)e 



2 



-1, 



(C4) 

(C5) 
(C6) 
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1 1 



Kgl(x,c) = JdnJ rfr2 e-( 1+c )l^-^l^( 1 + c )( 1 "^) + -( 1 -^) 







e -2(i+c)x + ( 3 + 8c + 4c 2) _ 2(1 + c)(2 + x + 2c(l + x))e~ 



41(i,-l/2)=A[i_( 1+l)e 



2(1 + c)(l + 2c)a; 2 



4 2 3 } (x,-l) = - 



1 -e" 



(C7) 



(C8) 



The coil functions involving three or four integrals are 

l l l 



K§1 (x, c) = J dn J dr 2 J dr 3 e x 



(1+c)|ti— T3|— x(l+c)|72— rs|+xc|n— T2| 



K§l(x, -1/2) = ^ x(l + e~ x ) - 2(1 - e"*) 





6 



and 



k£1(x,c 1 ,c 2 ) = J dn J dr 2 J dr 3 e 

ooo 
1111 

Kc\x,Ci,c 2 ) = dn dr 2 / dr 3 / dr 4 



£C3(1-T1 + |t 2 -T3|) + £C2(1-T2 + |ti-T 3 |) + :CC;l(1-T3 + |ti-T2|) 



(C9) 
(C10) 

(Cll) 



x e' 





xC3(\ti-T4\ + \t 2 -t :{ \)+xc 2 (\t 2 -T4\ + \t 1 -t 3 \)+xc 1 (\t 3 -T4\ + \t 1 -t 2 \) 



(C12) 



with c 3 = — 1 — C\ — c 2 . Clearly explicit expressions for these integrals can be obtained, these 
are already given elsewhere. 

The particular rod functions of interest are 



Si(y) 



K { n 2 \y,c) 

K m{y, c ) = I dsi j ds 
o o 



[-1 +cosy + ySi(y)], 



} } siny^/sf + 4 + 2sis 2 c 
dsi / ds 2 



o o 

i i 



yJs( + si + 2sis 2 c 



smyJsf + s|(2 + 2c) - Sis 2 (2 + 2c) 



W s i + s i(2 + 2c) - Sl s 2 (2 + 2c) 



(C13) 
(C14) 

(C15) 
(C16) 



and 



i i . / (:j) 

*3 



111 

K§l(y,c) = j dsx J ds 2 J ds 3 



ooo 

-(3) _ ,„ . 2 



(C17) 



T R1 = ( S l - S 3) + 2(si - S 3 )(s 2 - S3)C+ (s 2 - S 3 ) , 
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Ill 





i/v r S } 



T 



(3) 
R2 





,2 i /„ „ \2 , /„ „ \2 



S 3 + (Si - S 3 ) + (s 2 - S3) - 2s 3 (si - s 3 )ci 
2s 3 (s 2 - s 3 )c 2 - 2(si - S 3 )(s 2 - s 3 )(l + Ci + c 2 ), 



(C18) 



1111 

K^\y : ci,c 2 ) = J dsx j ds 2 J ds 3 J ds 4 



smyv4 4) 



r (4) 





\2 1 /„ „ \2 1 /„ „ \2 



= (si-s 4 ) +(s 2 -S4) +(s3-s 4 ) + 2(si - s 4 )(s 2 - s 4 )ci 
+ 2(si - S 4 )(s 3 - S4)c 2 - 2(s 2 - s 4 )(s 3 - s 4 )(l + Ci + c 2 ). 



Also, there are the so called mixed terms 

,(2), , 1 



4y — y cos y — 3 sin y — y Si(y) 
1 



(C19) 



(C20) 



= / <fei / ^ 2 



s 2 (2s 2 - si) 

r C5 1/ 



(Ac? - 3) sin(yr ( ^) + 3yr^ cos (yr^) 



.(2) 



(2). 



-(2h 



r (2) 

T cs 



S\ + S 2 — SiS 2 . 



(C21) 



Finally, we list the K function which appear as factors in vertices with nematic order 
parameters, 



JO S z 



1 d \ 



smys 



K 



(i) 

51 



1 

1 ds / 1 d 2 



dy 2 y dy / ys 
3y cos y — 3 sin y + y 2 Si(y) 
1 d \ sin ys 



(C22) 



y 2 y3 ^y ] yg 



8y5 



2/ 4 Si(y) + y(y 2 + 6) cosy + (y 2 - 6) siny 



(C23) 



*g/(v) 



1 ds f 1 d 3 Id 2 1 d ^ sin ys 



s 4 \ydy 3 y 2 dy 2 y 3 dy / ys 



1 



y 4 Si(y) + y(y 2 - 18) cosy + (18 - 7y 2 ) siny 



(C24) 



1 ds / d 4 6 d 3 15 d 2 15 d \ sin ys 



ys 



s 4 \ dy 4 y dy 3 y 2 dy 2 y dy 
^ [3y 3 Si(y) - 5y(y 2 - 42) cosy + 15(5y 2 - 14) siny 



and 



8y 



1 f 1 dsds' fid 2 1 d\ sin y(s-s') 



(C25) 



1 4 \y 2 dy 2 y 3 dy I y(s — s') 
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12y 5 



3y 4 Si(y) - 8y 3 + 3y(y 2 - 2) cosy + 3(y 2 + 2) siny 



/•i /-i rfsrfs' /id 3 1 d 2 1 d \ sin y(s-s') 
o io (s — s') 4 \y dy 3 y 2 dy 2 y 3 dy J y(s — s') 



% 5 



I/ 4 Si(y) + y (y 2 + 6) cos y + {y 2 - 6) sin y 



(2) /-i r±Jsd^_ (d^_Qd^ 

KsAy) -Jo Jo (s-s'Y Uv 4 V*/ 3 



6 <i 3 15 gP 15 d \ sin y(s — s') 



12y E 



) 4 \<iy 4 y^y 3 y 2 dy 2 y dy J y(s — s') 
9y 4 Si(y) - 64y 3 + y(9y 2 - 210) cosy - 15(y 2 - 14) siny 



and 



K s\y) = I dsi I ds 2 I ds 3 



(gl ~ g3)(g2 ~ 53) 

4 3)5 y 3 



(C26) 



(C27) 



(C28) 



T (3) 



S 2 + S 2 . + S 2 — S1S2 — ~~ S2S3- 



(l/ 2 4 3)2 - 3) sin(y4 3) ) + 3y4 3) cos(yr^ 3) ) 

(C29) 
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FIGURES 




FIG. 1. Parametrization of the configuration of a rod-coil diblock. 
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FIG. 2. The numerical value of the wave vector q* for TV = 40 as minimum of 
characteristic rod wave vector q* = 2ir /(Nfji). 
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FIG. 3. Phase diagram with ./V = 40, r = oj/x = 0. The critical point lies in the rod-rich region 
near f ro d ~ 0.73 and there is a triple point near f ro d ~ 0.42, separating the isotropic/BCC/HEX 
phases. 
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FIG. 4. Phase diagram with N = 40, r = uj/x = 1. The critical point still lies near f ro( i ~ 0.73. 
The three other dots are triple points. The dashed curves in the smectic C phase are contour lines 
for the angle 6, separating intervals of 5°. The smallest value for 9 is about 30° just at the phase 
boundaries near the triple points. 
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N=40, co/x=2 
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FIG. 5. Phase diagram with N = 40, r = u) j\ = 1. There is no critical point. All dots are 
triple points. The dashed curves in the smectic C phase are contour lines for the angle 9, separating 
intervals of 5°. 
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FIG. 6. Phase diagram with TV = 40, r = u/% = 4. All four dots represent triple points. The 
dashed curves and the thin solid curve in the smectic C phase are contour lines for the angle 0, 
separating intervals of 5°. 
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FIG. 7. Phase diagram with N = 40, r = lj/x = 8. There are no triple points. 
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FIG. 8. Phase diagram with N = r = w/x = 2. In this particular case there is no critical 
point and the dots are triple points. The dashed curves in the smectic C phase are contour lines 
for the angle 6, separating intervals of 5°. 
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